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Abstract 

Alternative gravitational theories described by Lagrangians depending on general functions of the Ricci 
scalar have been proven to give coherent theoretical models to describe the experimental evidence of the 
acceleration of universe at present time. In this paper we proceed further in this analysis of cosmological 
applications of alternative gravitational theories depending on (other) curvature invariants. We introduce 
Ricci squared Lagrangians in minimal interaction with matter (perfect fluid); we find modified Einstein 
equations and consequently modified Friedmann equations in the Palatini formalism. It is striking that 
both Ricci scalar and Ricci squared theories are described in the same mathematical framework and both 
the generalized Einstein equations and generalized Friedmann equations have the same structure. In the 
framework of the cosmological principle, without the introduction of exotic forms of dark energy, we thus 
obtain modified equations providing values of w e ft < — 1 in accordance with the experimental data. The 
spacetime bi-metric structure plays a fundamental role in the physical interpretation of results and gives 
them a clear and very rich geometrical interpretation. 
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I. INTRODUCTION 



In this paper we try to better understand and to analyze alternative theories of gravity 
depending on higher-order terms in the curvature invariant RU*") Rr^ v ), in relation with some 
very interesting and possible cosmological application and, in particular, in relation with their 
capability to explain the cosmological acceleration of the universe, both in early times (inflation) 
and in present time universes. Nevertheless we will focus our attention on the possible theoretical 
explanation of the present cosmological acceleration. 

Recent astronomical observations have shown that the universe is accelerating at present time 
< s ee Q md Q fa _ observation re8 * ; see Q ta the ob— about the anisotrop, 
spectrum of the cosmic microwave background (CMBR); see jj| for the results about the power 
spectrum of large-scale structure). Physicists have thus to face the evidence of the acceleration 
of the universe and should give a coherent theoretical explanation to these experimental results: 
a problem which up to now seems to be still unsolved! General Relativity in interaction with 
a perfect-fluid like matter and the cosmological principle, providing the standard cosmological 
models, fail to give by their own a theoretical framework to explain the acceleration of the universe. 
We are thus forced to introduce some kind of dark matter or dark energy, which are responsible 
for the acceleration of the universe, or to modify General Relativity such that acceleration is 
predicted (see for example |J|). 

Dark matter or dark energy models have been deeply investigated in relation with their capability 
of explaining the acceleration of the universe (see |(j and references therein), however up to now 
there are no satisfactory experimental evidences of the presence of the predicted amount of dark 
energy in the universe. The real nature of dark energy, which is required by General Relativity 
in this cosmological context, is unknown but it is fairly well accepted that dark energy should 
behave like a fluid with a large negative pressure. The dark energy models with effective equation 
of state w e ff (which determines the relation between pressure p and density of matter p) smaller 
than w e ff < — 1 are currently preferrable, owing to the experimental results of Q]. 
On the other side the simplest way of obtaining accelerated expansions within General Relativity 
is to introduce a positive cosmological constant Ml, an introduction which leads however to some 

n q 

theoretical and experimental problems and contradictions (see for example |5( and [7]]). We just 
want to stress here that models with a constant cosmological constant are not able to explain the 
evolution between different epochs of the universe, characterized by different values of acceleration 
(deceleration) . 
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The other possibility is to assume that we do not yet understand gravity at large scales, which 
means that General Relativity should be modified or replaced by alternativeeravitational theories 
of Gravity when the curvature of spacetime is small (see for example j|, U and references 
therein), providing modified Friedmann equations. Hints in this direction are suggested moreover 
from the quantization on curved spacetimes, when interactions among the quantum fields and the 
background geometry or the self interaction of the gravitational field are considered. It follows 
that the standard Hilbert-Einstein Lagrangian has to be suitably modified by means of corrective 
terms, which are essential in order to remove divergences [8j|. These corrective terms result to be 
higher-order terms in the curvature invariants such as R, R^R^, R^ ua ^ R^ va p, RO l R, or non 
minimally coupled terms between scalar fields and the gravitational field. It is moreover interesting 
that such corrective terms to the standard Hilbert-Einstein Lagrangian can be predicted in higher 
dimensions by some time-dependent compactification in string/M-theory (see [3]) and corrective 
terms of this type arise surely in brane- world models with large spatial extra dimensions 10]. As 
a matter of facts, if these brane models are the low energy limit of string theory, it is likely that 
the field equations include in particular the Gauss-Bonnet term, which in five dimensions is the 
only non-linear term in the curvature which yields second order field equations. In this framework 
Gauss-Bonnet corrections should be taken into account and cosmological models deriving from 

n 

the Gauss-Bonnet have been recently studied; see |12|] and references therein. 
As an alternative to extra dimensions, it is also possible to explain the modification to Friedmann 
equations (which could provide a theoretical explanation for the acceleration of the universe) by 
means of a modified theory of four dimensional gravity. The first attempts in this direction were 
performed by adding to the standard Hilbert-Einstein Lagrangian analytical terms in the Ricci 
scalar curvature invariant ^jj. A simple task to modify General Relativity, when the curvature 
is very small, is hence to add to the Lagrangian of the theory a piece which is proportional to 
the inverse of the scalar curvature or to replace the standard Hilbert-Einstein action by means 
of polynomial-like Lagrangians, containing both positive and negative powers of the Ricci scalar 
R and logarithmic-like terms. Such theories have been analyzed and studied both in the metric 
Q a» d ft. Pacini f „ rm a lismS Q, Q. It ^ th a t b o th b th e metr , a„ d tb. Pacini 
formalism they provide a possible theoretical explanation to the present time acceleration of the 
universe. Moreover a mechanism ruling the present dark energy dominance 1 (due to the universe 



1 Taking into account the transition of the universe from a decelerated era to an accelerated era, a scenario with 
w e ff transitting from values below —1 to values above —1 is actually preferable [lrj. 
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expansion) and the present cosmic acceleration has been proposed in this framework; see 
A discussion is open on the physical reliability of the Palatini and/or the metric formalism [lsl ] 
and on the physical relevant frame both in the metric and the Palatini formalism 19]. Up to now 
it appears that the aforementioned metric approach 2 leads to results which are in contrast with 
the solar system experiments and also that the relevant fourth order field equations suffer serious 
instability problems ^j]. On the contrary the Palatini formalism produces second order field 
equations which are not afflicted by instabi lity problems and are in acceptable accordance with 
the results of the solar system experiments jl5|. A discussion is actually open on the accordance 
of the Palatini formalism with the electron-electron scattering experiments 



The importance of modified theories of gravity depending on general analytical functions of 
the Ricci scalar is also related with the possibility of avoiding singularities in these cosmological 
models j^J and in the interpretation of black holes entropy in this context j^]. Recently, an 
explanation of the present day acceleration of the universe has been moreover formulated in the 



framework of non-symmetric gravitational theories 



determinant of the Ricci tensor 



23]. 



22 ] and in modified theories depending on the 



Encouraged by recent developments of cosmological applications of alternative theories of 
Gravity we consider in this paper Ricci squared Lagrangians in minimal interaction with matter, 



which have been deeply analyzed in 2J] in the vacuum case. As we already said before such 
Lagrangians are deeply related with quantum field theory: to remove divergences one has to add 
counterterms to the Lagrangian which depend not only on the Ricci scalar but also on the Ricci 
and the Riemann tensors 8]. It was proven in j^J that Ricci squared Lagrangians provide second 
order field equations in the Palatini formalism, such that the universality of Einstein equations 
and the universality of the Komar energy-momentum complex hold in vacuum. These remarkable 
results have important implications also in cosmological models. They imply that, in some 
sense, field equations for Ricci squared Lagrangians reproduce (apart conformal transformations) 
the standard Einstein field equations in the vacuum case, while in the presence of matter this 
equivalence might be broken. The geometrical structure of the spacetime manifold is very rich and 
it is endowed with an anti-Kahlerian structure, deriving directly from the variational principle of 
Ricci squared Lagrangians (see (3]). Spacetime turns out to have a bi-metric structure, or better 



We remark that field equations are in that case fourth order field equations. 
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a so-called metric compatible almost-product as well as an almost- complex structure with a Norden 
metric. The geometrical structure of spacetime is moreover characterized by a scalar-valued 
structural equation, which is simply obtained by contracting field equations with the metric j^J 
and controls the solutions of field equations. Lagrangians based on higher-order Ricci scalars 
which led to higher-order metric compatible polynomial structures have been considered both in 
purely metric and Palatini formalism in 



It was moreover shown in 



25l ] and that Ricci squared theories in vacuum give field equations 



equivalent to Einstein field equations with a cosmological constant, the value of which is fixed once 
the structural equations are solved and one particular solution of the structural equations is chosen. 
This fact is no longer true in the case of interaction with generic matter, where the solutions of 
the structural equation are dynamical (the same happens in the case of non-linear Lagrangians in 
the Ricci scalar; see |l4 | and references therein). The equivalence with Einstein field equations 
is hence broken and we obtain modified field equations, depending on the stress-energy tensor of 
matter involved in the theory. Nevertheless, we remark that field equations are once again second 
order field equations in the metric field. 

For cosmological applications we consider the physical metric g to be a Robertson- Walker metric 
and the stress-energy tensor of matter to be a perfect fluid one. In this particular framework 
deriving from the cosmological principle we obtain that the Levi-Civita metric h is conformal to 
the physical metric g, apart for a rescaling factor of the cosmological time. From our construction 
it follows however that the signature of h can be arbitrarily chosen (it can be either Riemannian, 
or Lorentzian or Kleinian, apart from some restrictions deriving from field equations). We are 
consequently able to introduce a generalized Hubble constant and modified Friedmann equations. 
A comparison with the f(R) theories is immediate. It is striking to notice that modified Friedmann 
equations are once more first order field equations, which prevent the appearing of instabilities 
as it has already been shown in the case of Ricci scalar theories in ^^|. This is an important 
consideration giving the Palatini formalism a deeper physical significance, in view of cosmological 
applications. 

An explicit example dealing with power Lagrangians in the Ricci squared invariant f(S) = (3S n is 
analyzed in detail and the Hubble constant is derived. It results that the deceleration parameter 
can be negative if particular values of n are chosen. Moreover we obtain values of w e ff that 
can be suitably fitted to the experimental results of [jj. Considerations are exposed about the 
frame changing, which means choosing h to be a FRW metric instead of g. Field equations and 
cosmological parameters are obtained and discussed also in that alternative (Jordan) frame. 
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The paper is organized as follows: we start in section 2 by considering the case of f(R) 
Lagrangians in a new matrix formalism, with the introduction of an operator P modifying the 
Einstein field equations We pass in section 3 to the more complicated case of Ricci squared 
f(S) Lagrangians and we analyze the field equations and the structure of spacetime in the 
case of interaction with matter. We proceed in section 4 with cosmological applications and we 
obtain modified Friedmann equations. In section 5 we discuss the relevant example of polynomial 
Lagrangians in the Ricci squared invariant. In section 6 we consider the theory in the alternative 
Jordan frame, where h is assumed to be a priori the FRW physical metric. 



II. COSMOLOGICAL MODELS IN f(R) GRAVITY 



We start considering non- 
and developed in [14| and in 



inear Lagrangians in the Ricci scalar invariant f(R), already treated 
in the vacuum case. We think that it is worth summarizing those 



theories in order to have a comparison with Ricci squared theories here developed and analyzed in 

n 

detail. We moreover modify the formalism introduced in ■14] to treat both Ricci scalar and Ricci 
squared theories in the same mathematical framework. 
The action for f(R) Gravity is introduced to be: 

A = A gTav + A mat = J ( v /det^/(i?) + 2 K L m ^))d 4 x (1) 

where R = R(g,F) = g a ^R a p{T) is the generalized Ricci scalar and R^(T) is the Ricci tensor of 
a torsionless connection T. The gravitational part of the Lagrangian is controlled by a given real 
analytic function of one real variable f(R), while ^fg denotes the scalar density | det || g^ v |||2 of 
weight 1. The total Lagrangian contains also a matter part L mat in minimal interaction with the 
gravitational field, depending on matter fields together with their first derivatives and equipped 
with a gravitational coupling constant k = 8irG. 

Equations of motion, ensuing from the first order a la Palatini formalism are (we assume the 
spacetime manifold to be a Lorentzian manifold M with dimM = 4; see 2sj]): 

f'(R)R iH (T) - \f{R)g^ = kT^ (2) 

V r a (^tetg~f(R)gn = (3) 

where T^ v = — -^= ^ mat denotes the matter source stress-energy tensor and V r means covariant 
derivative with respect to V. 



We shall use the standard notation denoting by R[^ v ) the symmetric part of R^ v , i.e. Ruv) = 
\{Rnv + Rv^))- I n order to get Q one has to additionally assume that L mat is functionally 
independent of T; however it may contain metric covariant derivatives V 9 of fields. This means 
that the matter stress-energy tensor T^ u = T^ u (g,^) depends on the metric g and some matter 
fields denoted here by ^, together with their derivatives. From (jHJ) one sees that y/detgf (R)g fMU is 
a symmetric twice contravariant tensor density of weight 1, so that if not degenerate one can use 
it to define a metric huv such that the following holds true 



v / drtg~f'(R)g^ = Vdrthh^ (4) 

This means that both metrics h and g are conformally equivalent. The corresponding conformal 
factor can be easily found to be f'(R) (in dimM = 4) and the conformal transformation results to 
be: 

Ku = f'(R)g^ (5) 

Therefore, as it is well known, equation @ implies that T = TLc(h) and R^ V ^(T) = R^ u (h) = R^. 
Let us now introduce a (l,l)-tensorfield P by 

P>f = g^R a u(h) (6) 

so that (j2j) re-writes as 

f{R)p; - \f{RW» = «zj (7) 

where, with an abuse of notation, T = = g^ a T av and from © we obtain that R = trP. 
Equation ((7J) can be supplemented by the scalar-valued equation obtained by taking the trace of 
(J7J), (we define r = trT) 

f'(R)R-2f(R) = K g a/3 T a(3 = KT (8) 

which controls solutions of ((7J). We shall refer to this scalar- valued equation as the structural 
equation of spacetime. The structural equation ©, if explicitly solvable, provides an expression of 
R = F(t) and consequently both f(R) and f'(R) can be expressed in terms of r. More precisely, 
for any real solution R = F(t) of (JHJ) one has that the operator P can obtained from the matrix 
equation 0: 

p= f(F{r)) I+ k f 



Now we are in position to introduce the generalized Einstein equations under the form 



g^P? = R»v ih) (10) 

where is given by 0a n d Pu is obtained from the algebraic equations © and © (for a given 
g^v and T^); see also jlJ| and For the matter-free case we find that R = F(0) becomes a 

constant implying that the two metrics are proportional and the operator P is proportional to the 
Kronecker delta. Equation is hence nothing but Einstein equation for the metric g, almost 
independently on the choice of the function f(R), as already obtained in [28j |. Also the standard 
Einstein equation with a cosmological constant A can be recasted into the form (|1U[). It corresponds 
to the choice f(R) = R — A. These properties justify the name of generalized Einstein equation 
given to (|10j) . In the presence of matter equation IjlflJI expresses a deviation for the metric g to be an 
Einstein metric as it was discussed in [1J|. It can be otherwise interpreted as an Einstein equation 
with additional stress-energy contributions deriving from the modified gravitational Lagrangian 
or possibly as a modified theory of gravity with a time dependent cosmological constant. 



A. Cosmological applications of first-order non-linear gravity 

We give here a brief summary of the results obtained in ^j], where we refer the reader for 
further details. We assume g to be the Friedmann-Robertson- Walker (FRW) metric which (in 
spherical coordinates) takes the standard form: 

1 



g = -dt 2 + a 2 (t) 



1 — Kr 2 



dr l + r 2 [dO l + sin 2 (0)dip' 



(11) 



where a(t) is the so-called scale factor and K is the space curvature {K = 0, 1, —1). We further 
choose a perfect fluid stress-energy tensor for matter: 



T^u = {p + p)u^u u + pg^ u 



(12) 



where p is the pressure, p is the density of matter and is a co-moving fluid vector, which in a 
co-moving frame (n M = (1,0,0,0)) becomes simply: 

In n o o\ 


pa 2 (t)r 2 
\ pa 2 (t)r 2 sin 2 (0) J 



P o 

n pa 2 (t) 



(13) 
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The metric h turns out to be conformal to the FRW metric g by means of the conformal factor 
f'(R), which can be moreover expressed in terms of r by means of (JHJ) and finally as a function 
of time b(t) = /'(i2(r)), by an abuse of notations. From IjlflJI we can obtain an analogue of the 
Friedmann equation under the form 



H 



a b \2 
a + 26/ 



K 

3b 



P + 



f(r) + kt 



2k 



K 
~7i 



(14) 



which can be seen as a generalized definition of a modified Hubble constant H = ( , taking into 
account the presence of the conformal factor b(t) which enters into the definition of the conformal 

n 

metric h (see [14J for details). This equation reproduces, as expected, the standard Einstein 
equations in the case f(R) = R. 

Considering the particular example f(R) = (3R n we have obtained that the Hubble constant for 
the metric g can be locally calculated to be: 



where: 



H A = 

r(n, w) 
s(n, w) 



sr(n, w)a 



-3fw+l) 



s(n,w)—z 



(15) 



2n 

3(3w-l)[3w(n-l)+(n-3)] 

2n 



-k(3w— 1) 



0{2-n) 



3iu(n-l)+(n-3) 

are functions of the power n and of the equation state of matter, through w. We remark that 
e = signi? = 1 for odd values of n and, on the contrary, e = ±1 for even values of n; see ^| for 
details. The deceleration parameter can be obtained from the Hubble constant by means of the 
following relation: 

H(t) \ 



q{t) :- 



1 + 



H 2 (t) 



(16) 



and from 1)15(1 it turns out to be formally equal to: 
q(t,w,n) = -1 + 



3(l+w) , \ - 

2n tr(n,w)a 



aa+w} 



s(n, w)Ka 



er(n,w)a n — s(n,w)Ka~ 2 



(17) 



_ 9 3(l + m) 

It follows that when the a term dominates over a 2n the deceleration parameter results to 



be positive, i.e. q(t, w, n) — » 0. On the contrary, when the term a 



3(l+m) 



dominates over a 2 (or 



in the case K = corresponding to spatially flat spacetime) the deceleration parameter results to 
be: 



l{w,n) 



-1 + 



3(1 +w) 
2n 



(18) 



which is negative for n < or n > 3 ( 1 +" J ) > owing to the positivity of (1+w) for standard matter; 
see This implies that the accelerated behavior of the universe is predicted in a suitable limit. 
In particular it follows that super-acceleration (q < — 1) can be achieved only for n < 0. The 
effective w e jf can be obtained (as in [ll|) by means of simple calculations from (|15|) and (|18|). It 
results to be, for this theory: 

2 , , 1 (w + 1) 

weff = p(n, «,)-- = -1 + (19) 

We remark that the range of —1.45 < w e ff < —0.74 for dark energy, stated in can be easily 
recovered in this theory by choosing suitable and admissible values 3 of n. We refer to jlj] for 
physical considerations and for more detailed discussions and examples concerning polynomial-like 
Lagrangians in the generalized Ricci scalar. 



III. RICCI SQUARED LAGRANGIANS IN MINIMAL INTERACTION WITH MATTER 
FIELDS 

We consider now the action functional: 

A = A gTav + A mat = J (Vdit^ftS) + 2 K L m ^))d 4 x (20) 

where S = S (g,T) = g^ a R^ au ^{T)g l/ ^ R^^{T) and R^ U {T) is, as above, the Ricci tensor of a 
torsionless connection V (see discussion after formula (|T])). The gravitational part of the Lagrangian 
is controlled by a given real function of one real variable f(S); see j^J]. Under the same assumptions 
of [3] and in 4-dimensional spacetimes M (dimM = 4) equations of motion ensuing from the 



variational principle in the Palatini formalism are 



3: 



2f'(S)g^R {lia) (T)R (M (T) - \f{S)g^ = kT^ v (21) 



Vl(^tetg-f(S)g^R {aP) (T)g^) = (22) 
where = — ^ mat denotes again the matter source stress-energy tensor. The above system of 



equations splits, as before, into an algebraic part (|21|) and a differential one (122(1 for the unknown 
variables g (the metric) and T (the connection). 



Ht! 
eg 



3 As already explained in |14 the parameter n should not be an integer, it can be any real number satisfying some 



reliability conditions; see 14] for further discussions and details. 
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28( (see also QJ) let us 



.14]) 



Following the general strategy elaborated for the matter-free case 24] and 
notice that ^/detg f (S)g fia R( a /3)9 is a symmetric (2, 0)-rank tensor density of weight 1 which we 
additionally assume to be nondegenerate. This assumption entitles us to introduce a new metric 
hav by the following definition 

Vd^hh^ = s/d^f'(S)g^R {a0) (T)g^ (23) 

The metric h is hence called a Levi-Civita metric since the field equations (|22(l and consequently 
r is the Levi-Civita connection of it: T = T^cih). The Ricci tensor of h can be simply defined 
as i?( Ml/ )(r) = Rfj_ u (h) = Rfj_ u . It should be easily recognized that eq. (|2*2*)) defines only up to 
multiplicative constant. Therefore the metric h is not a good candidate for a physically meaningful 
object. 

The algebraic equation (|21|) can be easily converted into the matrix form 

p2 -lk) m, + W(S) f (24) 

by using the endomorphisms P and T (i.e. (1, l)-tensorfields) as defined before: 

P = If = g» a R {au) (25) 
T = Tjf = g» a T av 

and I = 5u denotes the identity endomorphism, i.e. a Kronecker delta in dimension 4. In matrix 
notation one can also write P = g~ 1 R and T = g~ 1 T. 

Equation (|2*1)) can be supplemented by the scalar-valued equation obtained by taking the trace of 
J2U) or of (El 4 

f(S)S - f(S) = ^T aP ^ \t (26) 

which governs solutions of the matrix equations (|21() and 1)24(1 and we will define it as the structural 
equation of spacetime under analysis. We remark that in the vacuum case, as much as in the 
particular case of radiating matter (r = 0), we have that (1261) gives constant solutions for the 
values of S, so that the universality property of Einstein field equations still holds 24]. In the 
more general case of interaction of the gravitational field with matter we are considering, we will 
have that solutions of (|26|) are no longer constants, but they are related with the values of r. This 
means that the solutions of (|26j) are dependent on the choice of the stress-energy tensor for matter 



We remark that in this context 5* = trP 
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(at least on the trace of the stress-energy tensor) and moreover these solutions are dynamical, since 
r is generally time-dependent. 

The structural equation ((26|) can be formally (and hopefully explicitly) solved expressing S = F(t). 

This allows to reinterpret both f(S) and f'{S) as functions of r in the expressions: 

( f(S)=f(F(r))=f(r) 

{ 27 

1 f'(S) = f(F(r)) = f{r) 

where, for convenience, we will use in the following the abuse of notation f(F(r)) = f(r) and 
f'(F(r)) = /'(t). For any real solution S = F(t) of ()26|) it is hence possible to compute the 
operator P by solving the matrix equation (|24|) : 

^-M) i+ ^k) f (28) 

by simply calculating a square root of the endomorphism on the right hand side (under the condition 
/'(t) 7^ 0). The P tensorfield results consequently to be a function of P = P(t), due to (|27|) . 
Owing to the cosmological principle it results that r and consequently the operator P will be 
simply functions of time, once the stress energy tensor for matter is chosen. 

We remark that the solution proposed above for the matrix equation (|28|) is just one of the solutions 
of (|28|) and precisely it represents the simplest diagonal solution in the set of all possible solutions 
of (|28|) . The definition of P given in (|25|) should satisfy some restrictions, deriving directly from 
field equations of the theory (|21f) . as much as h should satisfy them. These conditions can be read 
in a differential form in (|25[) . which is however unsolvable, or translated into an algebraic expression 
(|28ft . This equation thus selects operators which are meaningful in the theory we are constructing. 
Off-diagonal solutions can also be found (as much as in j^]), but in the 4-dimensional case under 
analysis they are very difficult to be explicitly calculated. For our purposes we restrict thence 
overslves to the diagonal solution; more complicated solutions of these equations, in relation with 
the geometrical structure of spacetime, will be possibly analyzed in forthcoming papers. 
On the other hand equation (|23|) tells us the the metric h is conformal to a symmetric bilinear 
form; i.e., in matrix notations: 

h ~ (g-'Rg- 1 )- 1 = P- I g (29) 



Now we are in position to calculate the conformal factor, which results to be SI = -^2iOi[/ / (S')] 1 



and we will have in matrix notations that h = QP l g. Owing to the equations (|23j) and |2f 
respectively, it is then possible to set: 

det R = det g ■ det P (30) 
/ / (5) 4 detii = det/i 
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If we consider together the above equations (|3U[) it results we see that the conformal factor can 
be calculated to be f2 = /'(S^Vdet P, where detP can be simply obtained from the solution of 
equation (|28j) . once the structural equations (|26|) are solved. At this point we stress again that 
the conformal factor 0, is defined only up to an irrelevant multiplicative constant which has no 
influence on the physically measurable quantities g and T. 

We are thus able to express the metric h in terms of the operator P and the physical metric g from 
J21, as: 

V = = / , (r)VdetP(r)^ Q (P" 1 )" (31) 

where we have stressed the dependence of h from r, which follows from 1)26(1. Once again to obtain 
this expression for /i„„ explicitly we should require that (|27|) can be solved analytically. Having 
finally calculated Ptf and det P from the algebraic equations (|2*H|) and (|3T)|) (for a given g^ and 
T^jy) the generalized Einstein equation ensuing from (|21|) take the simple form: 

iV (h) = P« 9iia (32) 

with h^y given by 1)31)1 and now the physical metric g does not need to be Einstein. This expression 
for the generalized Einstein equations is formally the same obtained for non-linear Lagrangians in 
the generalized Ricci scalar in 1)10(1 . Differences arise in the definition of the operator P (compare 
expressions ® and (|28|)) and the metric h (compare expressions (J2J) and <)31j) ). which in this last 
case results, in general, to be no longer conformal to g. For the same reasoning as before one 
should easily realize that for the matter- free case equation (|32() becomes just Einstein equation for 
the metric h, with a cosmological constant depending on the analytical form of f(S). We remark 
once again that in the vacuum case we have that P is proportional to the identity and solutions 
of 1)26)) are constants. In the case of interaction with matter both P and f{S) = /(r) depend 
on the stress-energy tensor of matter, i.e. they are both dynamical. We thus skip from a static 
model equivalent to a standard Einstein theory with cosmological constant to a more complicated 
dynamical model, which is no longer analogous to Einstein Gravity. 



IV. FRW COSMOLOGY IN RICCI SQUARED GRAVITY 



For cosmological applications (as already explained in section III Aj) one has first to choose the 

physical metric, which is assumed to be g for the moment, to be the Friedmann-Robertson- Walker 

(FRW) metric, which (in spherical coordinates) takes the standard form (jll[) . i.e.: 

1 



g = -dt 2 + ar(t) 



Ll - Kr 2 



dr 2 +r 2 (de 2 + sin 2 (6)d^ 



(33) 
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Another main ingredient of the cosmological model is to choose the perfect fluid stress-energy 
tensor for matter, introduced in (j!2j) and in a co-moving frame in (|13j) . From the conservation law 
of the energy-momentum VT^j, = the consequent continuity equation takes the form: 



p + 3H(p + p) = 



(34) 



where H = ^ is the Hubble constant. The above continuity equation imposes standard relations 
between the pressure p, the matter density p and the expansion factor a(t) namely: 



p = wp 



p = rja 



-3(1+™) 



(35) 



with a positive constant rj > 0. As it is well-known the particular values of the parameter w £ 
{—1,0, ^} will correspond to the vacuum, dust or radiation dominated universes. Exotic matters, 
which are up to now under investigation as possible models for dark energy, admit instead values 
of w < —1 which are supported actually by experimental data 3]. We remark that the above 
expressions (|33)l imply that both p and p depend just on time, while they do not depend on space 
coordinates as an immediate consequence of the cosmological principle. This implies that the 
variable t is an implicit function of the cosmic time t. In order to find its explicit dependence of 
time one has to solve the Friedmann equation. 

From (|33f) and ()13|) it follows that T results to be, using the definition (|25[): 

/ -p \ 
p 
p 
\ p J 

All diagonal solutions of (|28|) can be thus calculated, using expressions (|33[) and IfMfiJ) : 



(36) 



p',; 



1 f(r) + 2 K p ^. / f(r)-2 K p 



(37) 



where we have formally expressed S = F(t) from (|27|) . where r = 3p — p. We introduce moreover 
= ±1; p, = 0, .. . ,3, ensuing from the square root of the operator P 2 . Notice that all possible 
choices of give rise to all possible diagonal solutions of the matrix equation but still corresponding 
to the same solution S = F(t). This exhibits a phenomenon of signature change in f(S) theories 
(see below and j^]). Reality condition forces us to assume that all three terms 



/'(r)/0, f(r)-2 K p^0 and /(r) + 2k P + 



(38) 
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have to have at the same time the same (negative or positive) sign. In what following we denote 
e = signf(r) = sign (/(r) - 2np) = sign (/(r) + 2np). 

It is hence possible to calculate the Levi-Civita metric h, which from IJ23JI turns out to be: 

V = \e^e\f'(T)\[(f(r) + 2Kp){f{r)-2Kp)]l x (39) 




//(r) + 2kp eia 2 2 2 2 2 2 \ 
" e0 V /(r)-2^ 'T^^' £2r ° Sm {0) ) 

where we denoted e = eoei£2€3- Neglecting an irrelevent multiplicative constant factor (which can 
be in general complex or imaginary) the above expression can be suitably rewritten, for convenience, 
as: 

V = 6(r)Diag (-e c(r), - , e 2 rV, e 3 r 2 a 2 sin 2 (#)^ (40) 

where: 



&0") = vW)l[(/W+2«p)(/(r)-2Kp)]i 

(41) 



c(r) = Tfp^ 



r)— 2np 

and 6(t) results to be a generalized conformal factor 5 between the two metrics g and h, while 
c(t) describes a rescaling factor for the cosmological time t. We notice that both the generalized 
conformal factor b(t) and the rescaling factor c(t) are positive definite by definition. The change of 
signature is related with coefficients e = ±1 and the freedom in their choice produces a multiplying 
of the Friedmann-Robertson- Walker manifold, which could be related with quantum cosmology 
phenomena. 

From the above expression l|40|) it is possible to notice that some choices of the value of e^, which 
is up to now completely free, will change the signature of the metric h, so that a signature change 
process appears as much as in j^]. If we choose all = ±1 to be equal we will obtain again a 
Lorentzian metric, at most with a different convention in signs. If any other choice is performed 
we will possibly have different signatures for the metric (corresponding to Euclidean, Lorentzian 
or Kleinian signaturs) and the t coordinate may then loose its preferred physical significance. 
The Ricci tensor of the metric h can be calculated from the expression (|40j) : it results to be diagonal 
with the following components: 



^oo - ^ 



n ab he be „ ( b\ 'b 'a 
-2-7 + 2— + 7- + 2 t " 2 7-4- 
ab a c be \ b b a 



(42) 



5 It is evident from the above expression that the two metrics h and g are no more conformal as they were in the 
case of the f(R) Lagrangians, apart from the very particular case of c(t) = const. However a suitable redefinition 
of the cosmic time variable restores the conformal relation between h and g 
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R 



11 



eoei 



c 4(1 -Kr 2 ) 



10- 



^d b 
a b 



a c 



a c 



be 
be 



+ 



R 
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(2- 



a 2 r 2 



c \ ab 



R 



£3- 



arr 2 sin 2 ((9) 




2 £c 6c + /a 
a c o c V a 



a c 




6 d\ 8K 
+ 2- +4- + ei-=- 
o a / or 



+ 



a c 




+ ei- 



8if 



+ sin 2 (fl)e 3 



The r.h.s. of the generalized Einstein equations (|32|) is obtained from (|37|): 

'f(T)-2 K p eia 2 



1 /(r) + 2Kp 
P,9au = ^J m Drag 



£2« r • £3« ?" sm (c/j 



2„2 . 



£1 ~ £2 
ei£2 



(43) 



/(r) + 2«p' 1 -ifr 2 ' 

Comparing expressions (|42j) and (|43|) we obtain that we must impose that e± = €2, which derives 
from simple algebraic consistent conditions on the generalized Einstein equations (|32j) . It follows 
that, like in the standard cosmological models, we have only two relevant field equations, corre- 
sponding to the 00 and the ii components. The values of en and €3 are completely arbitrary. We 
remark however the the choice of the value of 63 does not affect field equations as it cancels from 
field equations, as we will see later. Field equations are fixed once we have chosen the values of en 
and ei which modify respectively the 00 and the ii component of field equations. 
To obtain modified Friedmann equations, we have to take into account the relevant generalized 
Einstein equations, which are for the 00 component of Q32JI : 

a b 



a b 



a c 



a c 



be 
be 



b , d 
+ 2- +4- 
b a 



2eo f{r) -2k P 
3 V /'W 



(44) 



and for the generic the ii component 



6. 




ac be fa 

2 +8 - 

ac be \ a 



b d \ en 8Kc 
+ 2- + 4- + — — =- 



b a J ei a' 

Subtracting the first equation (fill from the second equation we obtain that the second 

derivatives of the scale factor o and of the conformal factor b both disappear and we get the 



modified Friedmann equations in the form: 

2 



a b 
a + 2~b 



ce jf{T) + 2Kp 
4 \ 



'/(t)-2kp Kc 
~ 6061 rf" 



(46) 



f(r) 12 V f'(r) 

where the expression on the l.h.s. can be defined as a modified Hubble constant (which is moreover 
analog to (|14|) : see also ^J]), which rules the dynamical evolution of the universe: 



a b 
a + 2b 



(47) 



6 We stress again that with the assumptions ei = 62 each ii component provides the same field equation, as it should 
be expected. 
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The r.h.s. of the modified Friedmann equations for Ricci squared theories differs however from the 
r.h.s. of (|14|) for Ricci scalar theories, as it should be reasonably expected. The evolution of the 
model is just dependent on the evolution of the scale factor a(t) and of the modified conformal 
factor b{t) (i.e. on H), while derivatives of the factor c(t) disappear. This fact is strongly analogous 
with the case of Ricci scalar theories. 

We remark that, as already observed before, the expression (J4"6|) and field equations depend only 
on the values of eo and e\. The sign factor eo appears as a constant in front of the r.h.s. of modified 
Friedmann equations, which can be rewritten as 



H 2 = e 



/(r) + kt + 2kp Kc 
6V/ , (^)[/(^)-2 = ^l" ei ^ r 



(48) 



so that a suitable choice of eo allows the r.h.s. to be always positive as expected. In fact eo has 
to be chosen in accordance with the prescription eo = sign f [r){f (t) + kt + 2np), provided the 
conditions are satisfied. However, we see that on the other hand ei appears only in the term 
related with the curvature K of the spacelike hypersurface. As it is obvious also from the explicit 
expression of h ()4U|) . choosing different values of ei is equivalent to change the sign of the spatial 
curvature. We remark finally that the choice of 63 is irrelevant for field equations. 



V. POLYNOMIAL LAGRANGIANS IN THE RICCI SQUARED INVARIANT 

We choose, as a relevant example to deal with, polynomial Lagrangians in S. In strict analogy 
with what has already been done for the Ricci scalar theories, polynomial Lagrangians can be 
considered as approximations 7 of any analytical expression in S in the suitable limit ^j]. It is 
hence worth investigating the behavior of cosmological solutions of Ricci squared theories described 
by means of Lagrangians which are pure powers of S: 

f(S) = (3S n (49) 

As a matter of facts polynomial Lagrangians can be approximated to pure power Lagrangians if 
the asymptotical behavior is considered and just the first leading term is taken into account. 
From the structural equations (|26|) we obtain that for the above pure power Lagrangian (|49|) in 

7 We are particularly interested in the cases of very small and large values of S, reproducing the cases of large and 
small curvatures of the universe, owing to the (linear) quadratic relation between S and R^ v . 
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the Ricci squared invariant: 



•^( T ) ~~ 2(n-l)' 

f'(r)- 



(50) 

2 (n-1) ' 

where we have to require n / 1 to avoid singularities in the theory (this implies that the case 
f(S) ~ 5 is not allowed 8 ). Since in the physically interesting cases one has r < we see that for 
generic n we have to assume 



P(n - 1) < 



Q 



However, for odd n we can allow (3{n — 1) > (see also |14( | in this context). Taking into account 
the standard relations: 



r = (3u> — \)p , p = wp , p = rja 



-3(1+™) 



(51) 



and performing straightforward calculations we obtain from (|41[) the generalized conformal factor 
(up to multiplicative constant): 



-3(l+™)(i-i) 



b(t) ~ a 



Performing further calculations by means of (|47|) it is simple to obtain: 



(52) 



H" 



(3w + l)(2n - 1) - 2 
4n 



n 2 



(53) 



We remark that in the particular case of w = —1 the expression (|53|) implies that H 2 = H 2 , 
independently on the value of n. Using the same relations, the r.h.s. of equation (|46[) results to be: 



H 



2 _ e (3w + l)(2n-l)-2 
6 y/n(3w- l)(3w + 3-4n) 



k(3w — 1) 
2(3(n - 1) 



j_ eo K Anw — (uu + 1) 



P 2 " --75 



ei a 2 y 3(1 + u;) — 4n 



(54) 



We stress that the rescaling factor c(t) is, in this particular example, independent on time: 

c(f) = 



4tun — io — 1 



(55) 



3w + 3 — 4n 

Combining equations (|53|) and (|54|) together we obtain that the Hubble H 2 constant for the physical 
metric g is: 



H 



3(l + m) gQ _„ 

eoP(n,w)a 2^ Q(n,w)Ka~ 

ei 



(56) 



This is similar to the presence of a methodological singularity for f(R) — R 2 in non-linear theories of Gravity 
depending on the Ricci scalar; see e.g. and 
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where we have defined: 

P(n, w) 
Q(n,w) 



8n 2 



3 v /n(3iu-l)(3iu+3-4n)[(3to+l)(2n-l)-2] 

2 



w(3tc-l) 
2/3(n-l) '/ 



i 

2 ii 



Anw— (w+1) 
3(l+«i)-4n 



4n 



(57) 



(3u>+l)(2n-l)-2 

From the above expression the deceleration parameter can be calculated by means of the standard 
formula already introduced in (j!6|) and it can be formally calculated from (|56j) under the form 9 : 

q(t, w, n) = -1 + m — ] — (58) 



J_ _3(1+™) 



Q(n, w)Ka 



-2 



We obtain consequently that when the a 2 term dominates over a" 
results to be positive, i.e. q(t,w,n) — * + , while when the term a 



3(l+w) 



the deceleration parameter 



30+w) 



dominates over a or in 



the physically very important case K = 0, the deceleration parameter will be: 

3(1+™) 



q(w,n) 



-1 + 



An 



(59) 



This implies that q(w,n) is negative for n < or n > 



3(l+w) 



> 0, owing to the positivity of 



(1 + w) > for standard matter. Taking also into account restriction on paprameters (w, n) 
coming from (|38|) the whole situation can be visualized on the phase diagram FIG. ^ 

Comparin g ex pression (|56|) and the standard relation which derives from General Relativity 



(see 



11] and 



it is easy to obtain that the effective value of w e ff deriving from Ricci squared 
alternative theories of gravity is: 

(l + w) 



w 



eff 



w,n 



-1 + 



2n 



(60) 



Both limiting values for w e ff(^,n) and w e ff(0,n) are marked on the phase diagram FIG.^ 



We can compare the values of q(w,n) and w e ff(w,n) obtained in this case with the values of 
q(w',n') Q18JI and w e ff(w',n') ()19j) obtained for alternative theories of Gravity with pure power 
Lagrangians of the Ricci scalar F(R) = f3' R n ' , already treated in 14]. It turns out that they differ 
just for a factor if the value of w is assumed to be fixed, which is equivalent to state that we are 
dealing with the same kind of matter. It is simple to see that: 

f q(w, n) = q(w',n') ^' 
{ w e ff(w,n) = w e ff(w',n') 
These results generalize and confirm the results already obtained in [31|] in the very particular case 
of quadratic Lagrangians. 



n 



(61) 



We say that the deceleration parameter can be formally calculated, as we do not know a priori if any physical 
solution exists in all cases considered. 
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FIG. 1: Phase portrait for the plane (w, n), where w is the fluid parameter characterizing matter and n the exponent 
of S in the Lagrangian. The shadow areas represent physically and mathematically admissible pairs. We notice the 
following: 

- for n > 1 we cannot have radiation (w = |) but dust is allowed; 

-for dust matter (w = 0) and n > 1, w e ff — > — 1 + , i.e., w e ff can be approach only from above; 

- in the contrary negative powers (n < 0) do not allow dust, although any dust-like matter can be allowed for large 
enough \n\. Moreover in this case w e ff — ► — 1~ is possible from below (super-acceleration). 

A. Polynomial Lagrangians 

As we already stated before, pure-powers Lagrangians in the Ricci squared invariant can be con- 
sidered as approximations of more physical polynomial-like Lagrangians of the type: 

= s+ n /wn + T^—S m (62) 

(1 + n)b n 1 — m 

(here both n > and m > 0, with 

We just consider for simplicity the case of flat universe K = 0. In the limit of small or large 
curvatures, corresponding to the cases of present time universe and early time universe, we obtain 
from the structural equations that the leading terms are respectively: 

f S — ► =^ -jj^- = KT 

\ S -> oo -f3S m = KT 
From (|59|l we deduce that polynomial Lagrangians provide an explanation for early time inflation 
assuming that m > | and they can provide an explanation to present time cosmic acceleration 
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assuming that some inverse power of the generalized Ricci squared curvature invariant is also 
present in the Lagrangian (i.e. a^O). 

This result generalizes previous results which have been obtained for Ricci scalar alternative theories 
of gravity Q|, Q], and they are related to the so-called Starobinsky inflation Q]. 

VI. CHANGING FRAME 

We have developed up to now a first order a la Palatini theory which after appropriate reduction 
turns out to be based, as we remarked before, on a bi-metric spacetime with an almost complex 
structure. In Section Hvl we have assumed g to be the "physical" FRW metric. However we do not 
know a priori which is the most appropriate frame in the bi-metric structured spacetime we have 
constructed in Section II I II This is the same problem already studied and examined in ^| , ^| 
and 19] in the case of Lagrangians depending on the Ricci scalar, where different frames result 
to be somehow inequivalent. We shall not comment here on this equivalence problem and refer 



the reader to a recent interesting discussion by Flanagan (|32| and ref.s quoted therein). In our 
understanding this important problem should be analyzed in more detail, also in relation with 
the physical consistency requirements and the presence of instabilities; we plan to treat it in a 
forthcoming paper (3!]]. 

In this framework it is thus worth considering also the case when the "Jordan" h is chosen to be 
the physical FRW metric, as we have already done in Q|. More precisely, according to (|iU)) we 
set: 

h = -e dP + ^ 2 (t)[ i ^ Kr2 dr 2 +r 2 (t l d6 2 + e 3 sin 2 (fl)V)] (63) 

to be, modulo signature, FRW metric with a new cosmic time t and a new scale factor A. The 
generalized Einstein equation (|l(Jj) can be also calculated in (t,x 1 )- coordinates. This is equivalent 
to the assumption that the metric h is the physical one (i.e., that we can use conformal Jordan frame 
instead of the original Einstein frame). In this case, one has to restore the standard Lorentzian 
signature by setting eo = t% = ±1. We consequently obtain: 

34 = ^^^) I = ^(4^) ! = ^m-* m 



A c(t)6(t) V 4/'(r) J 2|/'(t)| V /M + 2np J 2|/'(r)| 
for the 00 component while for the 11 component we find 

A ,A\* K _ 6q ( f( T ) + 2K P *y_ e / /(r) + 2/qA *_ e 1 
A + \A) + V b(r)\ 4/'(r) ) ~ 2|/'(r)| \f(r) - 2np ) ~ 2\f'(r)\ ° [T) 
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where A denotes now the differentiation with respect to the new cosmic time t. We have also taken 
into account that in this case df 2 = c{T)b{r)dt 2 and A 2 = b(r)a 2 . 



Now the analogue of the Friedmann equation takes the form 



H 



i 

e c2 



K 
A 2 



(66) 



with H 



12|/'(r)| * 

being the Hubble constant of the conformal metric h. Thus up to now arbitrary sign 



factor eo = ±1 can be adjusted as eo = sign (3 



in order to preserve the positivity. 



We specialize now to the case of pure power Lagrangians in the Ricci squared curvature invariant 
and in the meanwhile, as already stated before, to the case of polynomial Lagrangians in some 
suitable limit. Choosing, as already done in (|49|) . the Lagrangian to be: 



f(S) = (3S n 

we obtain respectively /(r) = 2 ( n -i) T anc ^ f'( T ) = n $ 
that the modified Friedmann equation in this frame is: 



H 



P(w,n) A 



A 



2/3(n-l) 

K 

" A 2 



(67) 

It follows from equation (|66|) 



where for convenience sake we have defined the coefficient P{rj,w,n) as: 



P(w, n) 

and the exponent A as 



|2n(l + 3w) -3(w + l) 



6/3 |n| \4wn 



A 



w - l|4|3u; + 3 - 4n|4 
12(w + l)(n - 1) 



2[3{n - 1) 

KT](3w — 1) 



(68) 



(69) 



(70) 



3(w + 1) - 2n(3w + 1) 
We can consequently obtain the deceleration parameter by means of formula (|58[): 

V ; PA X - KA- 2 V ' 

This implies that, in the limit when the term A~ 2 is dominating over A x we will have q(t) — > + . 
Otherwise in the limit when the term A x is dominating over A~ 2 or in the case K = 0, we will 
have 

An - 2>{w + 1) 



Thus 



which for big n behaves as 



q(t) -> q(w,n) 



w eff (w,n) 



An - 3(w + l)(2n - 1) 

2n(u> - 1) + w + 1 
4n- 3(u5 + l)(2n - 1) 



'u;, =too) 



1 — w 

3uT+T 
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VII. CONCLUSIONS AND PERSPECTIVES 



In this paper we have analyzed alternative theories of gravity depending on a Lagrangian assumed 
to be a general function of the generalized Ricci squared curvature invariant S constructed out 
of a dynamical metric g and a dynamical (torsionless) connection T. The Palatini formalism 
provides first order field equations for the metric and the connection T. A structural metric h 
is introduced, so that the connection turns out to be the Levi-Civita connection of h and h is 
consequently a Levi-Civita metric. A convenient spacetime bi-metric geometry is thus defined 
by means of generalized Einstein equations and it is controlled by means of structural equations; 
signature changing phenomena appear. This implies that the metric h can be either a Lorentzian, 
an Euclidean or a Norden metric, giving us an immediate and natural insight into quantum 
cosmology theories. 

To treat explicitly cosmological models we choose g to be a Robertson- Walker metric and the 
stress-energy tensor to be the stress-energy tensor of a perfect fluid. This allowed us to obtain 
modified Friedmann field equations and a modified Hubble constant related to a (generalized) 
conformal transformation factor b(t) between g and h and to a rescaling factor for the cosmological 
time c(t). The metric h can be considered to be FRW, too, so that it can be conveniently 
considered as a physical metric in place of the original g. Generalized Friedmann equations are 
obtained also in this framework. 

If we moreover specialize to the pure-power case f(S) = (3S n (with n an arbitrary real exponent) 
we have seen that, with suitable choices of the parameters involved, these models are able to 
explain the current acceleration of the universe. We obtain moreover that polynomial Lagrangians 
in the generalized Ricci squared invariant provide an explanation for the inflation of the universe 
in suitable limits 91. 



This paper was thus devoted to analyze the geometrical structure of spacetimes described by 
means of Ricci squared Lagrangians in interaction with matter; cosmological applications of this 
models have been analyzed, following the ideas of Q] and generalizing the effort to understand 



current acceleration of the universe in alternative theories of gravity and The relation 

between the geometrical bi-metric structure of spacetime (and in particular the signature change 
phenomena) and its cosmological implications is very rich in mathematical and physical significance 
and will form the subject of future investigations. 
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